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Fermion-number fractionalization without breaking of time-reversal symmetry was recently 
demonstrated for a field theory in (2 -I- 1) -dimensional space and time that describes the couplings 
between massive Dirac fermions, a complex-valued Higgs field carrying an axial gauge charge of 2, 
and a U{\) axial gauge field. Charge fractionalization occurs whenever the Higgs field either sup- 
ports vortices by itself, or when these vortices are accompanied by half-vortices in the axial gauge 
field. The fractional charge is computed by three different techniques. A formula for the fractional 
charge is given as a function of a parameter in the Dirac Hamiltonian that breaks the spectral 
energy-refiection symmetry. In the presence of a charge ±1 vortex in the Higgs field only, the frac- 
tional charge varies continuously and thus can take irrational values. The simultaneous presence of 
a half- vortex in the axial gauge field and a charge ±1 vortex in the Higgs field re-rationalizes the 
fractional charge to the value 1/2. 

I. INTRODUCTION 

The concept of fractional charge emerged from quantum field theory in 1976 when Jackiw and Rebbi showed that 
Bose fields can induce a fractional fermion number 1/2 for the relativistic fermions to which they couplcii The proper 
conditions for this mechanism of fractionalization are the following. First, the very notion of a fractional charge 
demands that the fermion number is a good quantum number. Second, the Bose fields must trigger the spontaneous 
breaking of a symmetry that opens up a gap in the singlc-particlc fcrmionic spectrum within the Born-Oppcnheimcr 
approximation. Third, the Bose fields must support local topological defects that nucleate single-particle fermionic 
bound states in the close vicinity to the defects. Fourth, this many-body quantum state is a finite energy eigenstate. 

The first requirement rules out mean-field descriptions of superconductors that can otherwise satisfy the remaining 
requirementsi^'^'^ The last requirement implies that the fractionalization of the fermionic charge is a long-distance 
and low-energy property of the many-body system, while the second and third ones insure a degree of robustness 
against local perturbations. This, in turn, suggests that the lessons learned from the quantum field theories in Ref. [l| 
could apply more generally to microscopic models encountered in solid state physics, thereby opening the possibility 
of a "table-top" measurement of the fractional charge. 

In fact, the work of Su, Schrieffer, and Heeger implies that the one-dimensional example of Ref. [l| can be thought of 
as an effective field theory that captures the relevant interactions between phonons and electrons in polyacetylenci^ 
Excitations with exotic quantum numbers (in relation to the fundamental electron constituents of the system), such 
as neutral objects carrying spin 1/2 or charge ±1 objects carrying zero spin, localize around a domain wall in the 
dimerization pattern of polyacetylene at the cost of a finite energy. Subsequent to this work, it was shown that 
exotic fermionic quantum numbers in one-dimension are not restricted to fractional values^i2ii2ii£ but can be tuned 
continuously by a small breaking of an energy-reflection symmetry assumed in Refs. [l| and[^, and defined below. 

With the discovery of the fractional quantum Hall effect, a different paradigm for charge fractionalization, one in 
which spontaneous symmetry breaking plays no role, was proposed by Laughlin for two-dimensional systems with 
strong breaking of time-reversal symmetryi ^^i^^ Central to this paradigm is the notion of topological order, a global 
property that characterizes an otherwise featureless incompressible liquid state of matter by the finite degeneracy of 
the ground state if the system is defined on a surface of non-trivial topology, with the degeneracy depending on the 
genus of the surfacCfi^ The fractional charge is intimately connected to this ground-state degeneracy, leaving no room 
for a continuoulsy varying fractional charge and, in particular, for an irrational charge. Since then, the preferred 
route towards charge fractionalization without time-reversal symmetry in two and higher dimensions has occulted 
any mechanism based on sponteneous symmetry breaking, presumably because it is believed that the energy cost for 
fractional charges is prohibitive in all but one dimension. 

However, as a matter of principle, this need not be so as was already shown by Jackiw and Rebbi in three- 
dimensional space when coupling Yang-Mills fields through the minimal coupling to Higgs fields and to Dirac fermions. 
Quantization of the Dirac fermions in the static background of a t'Hooft-Polyakov monopole nucleates a fcrmionic 
bound states with the fractional charge 1/2 at a finite cost in energy 

Of course, one might object that this three-dimensional example of charge fractionalization is unlikely to be realized 
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FIG. 1: (color online). The honeycomb lattice is shown in (a) and (b). The honeycomb lattice has 2 interpenetrating Bravais 
sublattices colored in blue and black, respectively. The electronic hopping amplitude is enhanced on thick red bonds while 
it is reduced on the thin yellow bonds relative to the magnitude t of the nearest-neighbor hopping amplitude. The so-called 
Kekule pattern dimerization pattern in (a) opens an energy gap for the single-particle fermionic levels and it maps to the 
complex- valued Higgs field of Ref. [S^ in the continuum limit. A dimerization pattern that shifts the relative separation of the 
Dirac points is shown is shown in (b). In the continuum limit, it maps to the axial vector potential introduced in Ref. [2ll. 

on the energy scale of the electron volt that governs solid state physics, a prerequisite for a table-top measurement 
of charge fractionalizatiouFii We do not know of a realistic three-dimensional model for band electrons coupled to 
bosonic collective modes that mimics Dirac fermions and Higgs fields coupled with each others and coupled minimally 
to Yang-Mills in the continuum limit. In two and three dimensions band-theory generically predicts an insulating or 
a metallic state of matter. In one-dimension the Fermi surface is generically realized by an even number of discrete 
points, thus providing the low-energy and long-wave-length limit of the tight-binding model with a Dirac structure 
for free. 

Semimetals, the most famous example of which is graphite, are exceptions to the hegemony of the band-insulating 
and of the metallic states of matter. Graphite is made of sheets of graphene, a honeycomb lattice made of carbon ions 
bound through sp^ orbitals, and where the fourth valence electron of each atom lazily revels predominantly between 
planar nearest-neighbor sites. The Fermi surface at half- filling for an isolated graphcne sheet is made of two isolated 
points.— The excitation spectrum around these two Fermi points endows the band electrons with a four-component 
Dirac structure owing to the Nielsen-Ninomiya theoremii^ 

Although this example of fermion-doubling is often viewed as a curse for the realization of quantum anomalieS j^^'^^i^^ 
it is this very property that opens the door to charge fractionalization without the breaking of time-reversal symmetry 
through spontaneous symmetry breaking, as shown by Hou, Chamon, and Mudryi^ The real-valued static fluctuations 
depicted in Fig.[lja) about the uniform nearest-neighbor hopping amplitudes of graphene are, in the continuum limit, 
represented by a complex- valued Higgs field that interacts with the four-component Dirac fermions. In the Born- 
Oppenhcimer approximation, a constant value of this complex-valued Higgs field breaks spontaneously an effective 
axial C/(l) symmetry of the continuum limit and opens up a gap in the single-particle fermion spectrum. If the phase 
of this complex- valued Higgs field is defective in that it carries a vortex, it nucleates single-particle mid-gap states that 
carry the fractional charge ±1/2 per state. The energy cost is not finite, however. In the continuum approximation, it 
grows logarithmically with the separation between the vortices. On the lattice, it even grows linearly with the vortex 
separation if the wave vector of the fluctuation of the hopping amplitude is commensurate with the reciprocal lattice. 

Jackiw and Pi showed that the energy cost of a vortex in the complex-valued Higgs field can be made finite if the 
complex-valued Higgs field and the Dirac fermions couple minimally to two real-valued Bose fields that realize the 
vector components of an axial gauge field and if this axial vector gauge field also supports vorticesi^ A honeycomb- 
lattice regularization of an axial gauge field without a vortex is shown in Fig. [IJb)i^ Alternative realizations of 
an axial vector potential also arise when the graphene sheet is curved^^i^ii^ or wrapped into fuUcrenesr^i^l into 
nanotubes^Sii^ and about a conei^i^ 

Charge fractionalization in one dimension can be tuned continuously by breaking an energy-reflection symmetry, 
which is defined below. This property survives in two dimensions and gives a mechanism for charge fractionalization 
that is fundamentally different (and thus potentially observable) from the mechanism for charge fractionalization 
that relies on topological order. By relaxing the condition of a finite energy to that of a logarithmically diverging 
energy, it was argued in Ref. [13 that a small staggered chemical potential that distinguishes carbon ions sitting on 
the now nonequivalent triangular sublattices of the honeycomb lattice can make the fractional charge irrational. This 
irrational charge was calculated analytically in the continuum limit and numerically for a lattice regularization in 
Ref. [2^. Remarkably, it was also found that the condition for finite energy in the continuum limit, i.e., the presence 
of a vortex in the axial vector potential, removed any dependence of the fractional charge on the staggered chemical 
potential. 

The purpose of this paper is to give three detailed and alternative derivations of the fractional charge that supplement 
the derivation from Rcf. |22|. Graphcne is of course not the only road to a semimetal in two dimensions. Threading the 
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elementary plaquettes of a square lattice with half a flux quantumai also realizes two nonequivalent Dirac points at 
half-filling. The Higgs field is then realized by a columnar pattern of dimerization whereas the axial vector gauge field 
is realized by a staggered pattern of dimerization.— It was shown numerically in Ref. [13 that the Z4 vortex defined 
by the four possible columnar patterns occupying the four quadrants of the square lattice pins the fractional charge 
±1/2 at the site where the four columnar patterns meet. The Z4 vortex is a discontinuous version of the vortices 
studied numerically in Ref. [l^. The fact that different lattice regularization of vortices carry the same fractional 
charge illustrates the fact that the fractional charge is independent of the short-distance regularization. This property 
will become obvious in the analytical calculations of the fractional charge that we are going to present. 

The paper is organized as follows. The quantum field theory is defined in Section [TTl The charge induced by the 
vortices in the Higgs or axial gauge fields is related to the spectral asymmetry in Section Hill The spectral asymmetry 
is computed in Setion IIVI The fractional charge is computed a second time by a perturbative expansion of the Dirac 
propagator in Section|Vl Finally, the fractional charge is computed on the basis of symmetry arguments in Section lVll 
A summary comprises the last Section IVlII 

II. DEFINITIONS 

In recent Letters j^i^i^^ we have developed the theory of charge fractionalization in planar models with topological 
defects encoded by vortices. The models are elaborations on graphene, with dynamics linearized around two Dirac 
points (the two inequivalent points in the first Brillouin zone at which the conduction and valence bonds of graphene 
meet). In a familiar fashion, the Schrodinger equation for the Bloch states at low energies and long wave length 
measured relative to the Dirac points takes a Dirac-like 4x4 matrix form, for a four- component "spinor" , which 
interacts with further scalar and gauge fields. The scalar field and the gauge fields are induced by fluctuations in the 
hopping amplitudes of the underlying microscopic tight-binding model. In this Section, we start with definitions. 

In second quantization, the planar Hamilton density reads 

n = i^^[a-{p- 75A5) + li - i75<^2) + R^\ V- = ^^H^. (2.1) 

Here, ?/;^ and 4' ^re creation and annihilation operators for four-components Dirac fermions, respectively, p = 
—i{d^,dy); is an axial vector gauge potential {Af,A^); (p^ and are the real and imaginary parts of a com- 
plex scalar field ip — (pi + 'np2 ; and /i is a field that acts like a staggered chemical potential - the staggering is governed 
by the matrix R. All fields depend on the three- vector = {t,r) = {t,x,y). The matrices in H are conventional 
4x4 Dirac matrices: 

« = (a^a^)^ (3^(^[iy 75 ^ = ( J i*! ) , (2.2a) 

where the "third" a-matrix. 




(2.2b) 



participates in the definition of 75 and also coincides with the matrix R = . (The matrices (^123 ^-re the standard 
Pauli matrices.) 

The Lagrange density corresponding to Eq. (|2.ip . 

£ = iij''dtij-H, (2.3a) 

is presented in covariant notation as 

C = ^Y (ia, + 75^5.) - (^1 - i75'^2) - 7V] V', (2.3b) 

with 

7° ^ /3, ^ ^ V^7°, 7 ^ /3a = ( ° "o*^ ) , ^'^m^0a^=(^ ) . (2.3c) 

[In Eq. (|2.ip . we set the axial scalar gauge potential Al? to zero.] 
Two gauge transformations leave the model unchanged. There is the local axial gauge symmetry 

^^e'-^^V, A,,^A,, + d,uj, ifi-^e^'^ip, ^i ^ ^x, (2.4) 
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where u; is a real- valued field and the index v = t^x, y. Also, there is the global phase symmetry, 



A 



(2.5) 



where nj is a real-valued number, which acts on the four components of the spinors. The latter leads to the conserved 
fermion (charge) number current. 



We shall show that the charge 



Q 



(Pr p{r) 



(2.6) 



(2.7) 



fractionalizes when the background bose fields are topologically nontrivial. 

The model possesses the usual discrete symmetries under the parity transformation V defined by 



f {t,x,y) 
ipit,x,y) 
A'Ent,x,y) 

^{t,x,y) 
I M(i,2;,?/) 



{t,-x,y), 

ij^ji'ip{t,-x,y), 

A5^{t,-x,y), 

~A-{t,-x,y), 

^{t,~x,y), 

-p{t,-x,y), 



(2.8) 



the charge conjugate transformation C defined by 



A - 

















(2.9) 



and the time-reversal transformation T defined by 



r : 



{t,x,y) 

il}{t,x,y) 

A'iit.x.y) 

v{t,x,y) 

p{t,x,y) 



{-t,x,y), 

j^'-f^ip'<{-t,x,y), 

A'^{-t,x,y), 

ip*{-t,x,y), 

fJ-{~t,x,y), 



(2.10) 



where one should remember that T is antiunitary so that complex conjugation of coefficients is implied. 
The theory possesses another discrete symmetry 



(2.11) 



In the lattice (the honeycomb lattice relevant to graphene, for example), the definition of parity depends on the axis 
used for the reflection; the transformation in Eq. (|2.8p corresponds to a reflection with respect to an axis that cuts 
through the bonds of the honeycomb lattice. 

When the staggered chemical potential p is dropped, i.e., the last term in the square brackets of Eq. (|2.ip or 
Eq. (|2.3bp is absent, the matrix R anticommutes with the remaining matrices in the single-particle Hamiltonian H of 
Eq. (|2.ip . Therefore, R maps positive energy eigenfunctions to negative energy eigenfunctions "i! _^ and vice- versa, 
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R^, 



(2.12) 



We call this an "energy-reflection symmetry" . 

We shall examine the Dirac theory with a specific vortex configuration for the Bose field (f, taken as a static back- 
ground, and with another specific vortex configuration for the axial gauge field Ag, also taken as a static background. 
The polar decomposition of the scalar field Lp is 



(r)e 



me 



= arctan 



y 



(2.13) 
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where the magnitude (f) oi ip vanishes at the origin (f>(r = 0) and tends to a nonvanishing ^(cxd) for large r. The integer 
n measures the vorticity encoded by the singular nature of the phase of the complex field tp at the origin. The axial 
gauge potential vanishes in the time component 

A°ir) = 0, (2.14a) 

while the spatial component reads 

Alir) ^ -ne/^a,{r), (2.14b) 

where a^{r) vanishes at the origin and tends to 1/2 at large r. The line integral over Eq. (|2.14bp along any closed 
curve that encircles once the origin yields the same number, proportional to the vorticity n. Finally, the chemical 
potential fi, also taken as a static background, is without topological structure and achieves a nonvanishing value 
/i(oo) at infinity. We shall take fj, to depend only on r, but it could also be constant. 

In the absence of the staggered chemical potential, the Dirac equation possesses \n\ zero-energy, normalizable 
solutions. These are the mid-gap states, eigenstates of R. Mostly, we consider the n = —1 case, with a single mid- 
gap state, ^fp, which remains bound even in the presence of the axial vector potential; turning on the axial vector 
potential changes the wave function profile, but the zero eigenvalue remains. We assume that there are no other 
bound states. When the staggered chemical potential is present, but never very large, the mid-gap state migrates to 
a shifted eigenvalue; however it still remains isolated in the gap. 



III. QUANTUM MECHANICAL ANALYSIS 



The following argument shows that without the staggered chemical potential fi the charge is —1/2 when there 
is a single normalizable mid-gap state '^q that is unoccupied. (When this mid-gap state is occupied, the charge is 
— 1/2+1 = +1/2.) The charge density arises from filling the negative energy continuum states of the Dirac equation. 



p{r) 



dE [^l{r)^E{r) - Tyr)T^(r 
dE hUr)'fEir) - TUr)T^(r) 



(3.1) 



where the second equality follows from the first due to the energy-reflection symmetry present in the problem at ^ = 0. 
The quantity T^T^ is constructed from reference states which solve a Dirac equation with a topologically trivial 
background and also possess the energy-reflection symmetry. In other words, the topologically determined charges 
that we compute are measured relative to a reference charge of a system with a topologically trivial background 
and possessing the energy-reflection symmetry. This procedure is needed to remove infinities. The reference wave 
functions form a complete set. The continuum wave functions ^ in the presence of the vortex - we call them 
the vortex states - are not complete; the mid-gap state is missing 



OO OO 

J dETUr)TE{r')^ j dE^l^r)^ eW) 



*J(r)vI'o(r'). 



(3.2) 



It therefore follows from combining Eq. p.ip with Eq. (|3.2p that 

p{r) = -^1'S(r)*oM 



and 



Q = j dVp(r) = -i. 



(3.3) 



(3.4) 



In the presence of the staggered chemical potential /x, the energy-reflection symmetry is no longer available to pass 
from the first to the second equality of Eq. (|3.1|) . However, we may proceed as follows. We suppose that the reference 
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states still possess the energy-reflection symmetry, so in Eq. (|3.ip we may still use this symmetry for thcmi^ 







J di?T^T^ = i J dETlTE = l J dE^i^E + l^l^b (3.5) 

— oo — oo — oo 

The last equality is again the statement of completeness of the continuum reference states and the continuum vortex 
states supplemented by the isolated bound state ^f^^, which is no longer at zero energy but has migrated to some other 
value in the gap between the continuum states. Using Eq. (|3.5[) in the first equality of Eq. (|3.ip leaves 

oo 

Pir) = -l^l{r)%ir)-^ J di? sign(i?) vl/t^(r)vl/^(r). (3.6) 

— oo 

It remains to evaluate the remaining integral, which is recognized as the "?7-invariant" , also called "spectral asymme- 
try". Note that with energy-reflection symmetry the integral vanishes, leading to the previous result (j3.3[) . In the 
above derivations, it is assumed that the "vacuum" is defined with the mid gap state unoccupied, and furthermore 
that the migrated state $j has positive energy so that it remains unoccupied in the definition of the vacuum. If the 
mid-gap state is occupied and/or the migrated state has negative energy, there occurs a sign change in Eq. (|3.3p that 
affects to the first term of Eq. (|3.6p . 



IV. SPECTRAL ASYMMETRY AND FRACTIONAL CHARGE 

We begin by putting the Hamiltonian Eq. (|2.ip in a more convenient form. This is done with the following unitary 
transformation 



H ^ H = THT-^, (4.1a) 



where 



The result for H is 



T ^ ( '^^^ ^+] , = Tt, -±^\{l± -3). (4.1b) 



^=(pt:) (4-2a) 



where the differential operator T) and its adjoint T>^ are given by 

V = ia'id. + e.^AD+iip^+a^if^, (4.2b) 

2?t = ia'(a,-e„-A^5)-i(/Pi+ (73(^2- (4.2c) 
With the factorization of the time dependence {u,v) — e^'^* {u^,v^), the stationary Dirac equation reads 



•nt M ) - ^ (4-3a) 
or in terms of components, 

Vve = {E + ^i)ue, (4.3b) 

VhiE^iE- fi)vE. (4.3c) 

Generally, fi can be a function on space-time. In the remainder of this Section, we shall set it to be a positive constant; 
the other background fields, ip{r) and Ai^{r) arc position dependent and static, with asymptotes quoted in Eqs. (|2.14p 
and ((TT31) . 
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A. Zero-mode solutions 



It is convenient to begin by considering two special cases, particular solutions of Eqs. (j4.3bp and (|4.3cP where either 
E ~ II OT E = —fi. These solutions would become zero modes of the Hamiltonian when /i = and they play a special 
role even when /i 7^ 0. 

Let us begin with the case where the energy eigenvalue E = ji. Then, from Eqs. (j4.3bp and (|4.3cp . it follows that 

V^u^{r) = 0, (4.4a) 

First we observe that if v^{r) were identically zero, u^{r) would also vanish and there is no solution. So we assume 
that v^{r) 7^ 0. Then, operating with on Eq. (j4.4b|) and using Eq. ((i:ia|) yields V'^Vv^{r) = 0. The following 
argument implies that Vv^ir) = 0. Consider 

= j St vl{r)V^Vv^ir) = J d\ \Vv^{r)\\ (4.5) 

Here, we are assuming that the spinor v^{r) obeys boundary conditions so that the differential operator is indeed 
the adjoint of I?, i.e., surface terms produced by partial integrations in the intermediate steps in (j4.5p vanish. Since 
the last integral vanishes, its positive semi-definite integrand must also vanish and wc conclude that 

Vv^Xr) - (4.6) 

while Eq. (|4.4ap implies that u^{r) = 0. 

Thus, we find that, when Ea. (|4.6p possesses a normalizable solution, there exists a positive energy bound state with 
E^li, 

pC) (.,>)) ^ "(.,>))• (") 

Similar reasoning establishes the occurrence of a negative energy bound state with E = —fi when there exists a 
normalizable solution of the equation 'D'^u_^(r) = 0, 



-,1 V\ (u_^{r)\ _ (u_^{r) 



dMu^,{r)? = 1. (4.8) 



The existence of solutions of the equations 'D'^u{r) = and T>v{r) = and the number of solutions of each kind 
are determined by the topological properties of the background fields, A^{r) and (p{r). An index theorem implies 

Index(i? ) = dim kcr V - dim kcr ^ n (4.9) 

where ker denotes kernel and n is the vorticity defined in Eq. (|2.13p . The implication of this index theorem was seen 
explicitly in Ref. d where solutions of 'D^u{r) = and I)v{r) = were constructed for the case of the highly symmetric 
profile of the vector and Higgs fields given in Eqs. (|2.14p and (|2.13p . It was found that, for a given vorticity, either 
one or the other of these equations has solutions, not both. Which equation had solutions depended on the sign of n. 
It was argued that the number of solutions is given by |n| and, when n = ±1, the solutions were found explicitly. 

The proof of the index theorem (|4.9p was given in Ref. l35l . The index theorem counts the difference indicated 
in (|4.9p . It proves that this is so, independent of the details of the profile of the vector and Higgs fields but with 
the assumption that, whatever they are, they are obtained by smooth deformations of the symmetric configurations 
in Eqs. (|2.13p and (|2.14p . Ref. |33 also presented a proof of a vanishing theorem, that either dim ker I? = or 
dimkerPt = 0. Combined with the index theorem, it implies that 

n>0 : dimkerPt^O, dimkerP^^n, (4.10a) 
n<0 : dim ker ^ 1 71 1, dimkerP = 0. (4.10b) 

A computation of the spectral asymmetry of the Hamiltonian in a spirit similar to the one that will be given in 
the remainder of this Section was originally presented in Ref. [s^. Equation (6.29) of that paper contains a result for 
the spectral asymmetry from which the index can be deduced by taking the parameter k (our fj,) to zero and which 
agrees with Eq. (|4.9p above. The general formula for the spectral asymmetry in their equation (6.29) also agrees with 
what we shall find in the following. 
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B. Non-zero mode spectrum 

Now, we shall look for eigenspinors of the Dirac Hamiltonian which do not have eigenvalues E = ±ii. 
From Eq. (j4.3cp we can solve for the lower components of the spinor in terms of the upper components 

v^ir) = -J—V^u^ir). (4.11) 

ii — /i 

Then, using (|4.3b[) we see that the upper components must obey the Schrodinger equation TyD^u^{r) = 
(i?^ — y?) u^{r). To find solutions, we begin with the eigenvalue problem 

W^u^ir) = \u^{r), A > 0. (4.12) 
We assume that we can find a complete orthornormal set of solutions of this equation, 

f d'ru{{r)uy{r)^5^y, J2u,{r)u{{r') ^ 6{r ~ r') 1. (4.13) 

Generally, the spectrum will contain both bound and continuum states. For continuous spectra, the right-hand side 
of the first equation above should be replaced with a Dirac delta function and the summation on the left of the second 
equation should be replaced by an integral. We shall assume that these replacements, where needed, are understood 
in Eq. (|4.13p . We can use the two-component spinor u-)^{r) to construct a normalized four-component spinor which 
solves the stationary Dirac equation. For each eigenvalue A, 

VaT7?-/A^/ u^{r) 



For every Uy^ which is a solution of the Schrodinger equation (j4.12p with positive eigenvalue A > 0, we obtain two 
solutions of the Dirac equation, one with positive energy, E = yX + Jj? and one with negative energy E = — •\/a + /?. 
Unlike the zero modes that we discussed in the previous Subsection, where there was either a positive or a negative 
energy solution, here, the positive and negative energy solutions of the Dirac equation are paired: for each positive 
energy solution there is a negative energy solution and vice- versa. This implies that, if there are bound states other 
than the zero modes, they must occur in positive and negative energy pairs. Thus, bound states, other than the zero 
modes, will not contribute to the spectral asymmetry. We will see this explicitly in the following. However, for states 
in the continuum spectrum, the pairing tells us only that the spectrum occurs symmetrically about zero: for example, 
there is continuum spectrum in the symmetrically placed regions E > ^w? -I- and E < — ^rn? + y? . It does not 
tell us about the density of states in these regions, which can still be asymmetric. 

C. Charge density 

Let us examine the charge density of the ground state of the system that we are considering. The charge density 
is given in Eq. p.6p 

= - -\Y^ sign(S)«'tj(r)1'B(r). (4.15) 

Here, we have included both types of zero modes. Depending on the sign of the vorticity, only one of them will 
be non-zero and will have mutliplicity given by the magnitude of the vorticity. A sum over these degenerate wave 
functions is implied in the first two terms on the right-hand side of (|4.15p . We have also assumed that [i is positive, 
so that V ^ is a positive energy state and is a negative energy state. We shall restore the possibility that /i could 
have a negative sign later, where it will simply lead to a flip in sign from the contribution of the zero modes. 
Now, using Eqs. (|4.14ap and (j4.14bp . we find that the third term in the right-hand side of (|4.15p is 



^— ^= f4(r)^,(r)-i(l?t^,(r))^2?t^,(r)) . (4.16) 



A>0 
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Using the fact that u-^^ satisfies the Schrodinger equation (|4.12p leads to 
p{r) = l:ul^{r)u_^^{r) - lvl{r)v^{r) 



2 ' ' 2 

-Y,—-jL== [ulir)VV^u,ir) - {V^u.ir))^ V^u.ir)) . (4.17) 



A>0 

The last terms in this expression are a total derivative 

p{r) = ^ul {r)u_ {r)^l-vl{r)v (r) 



2 "-/jv 

+d ■ ^^^l^^^, i'T -D^^xir)) ■ (4.18) 



A>0 



The total charge is a volume integral of the charge density. If we volume integrate the last term in the equation above 
and use Gauss' theorem, it will be expressed as a line integral on the circle at infinity of the quantity which is to 
the right of the derivative operator. Thus we see that the charge will depend on the asymptotic form of the wave- 
functions. We observe that, consistent with our discussion after Eqs. (|4.14ap and (j4.14bp . since the wave-functions of 
bound states fall off exponentially at large distances, bound states will not contribute to the charge. Only continuum 
states are important. Further, studying the asymptotics of the continuum states will allow us to compute the total 
charge. What will make the task easy is the fact that the volume integral of the part of the last term in (|4.18p will 
pick up contributions which go like 1/r. 

Before wc do that, we re-organize the expression for the charge density. We use the identity 



/Lt did fl 



and the Schrodinger equation (|4.12p to re-write Eq. (|4.17p as 
p{r) = ^u^^(r)u_^(r) - ]^vl{r)v ^,{r) 



2 



(4.19) 




(4.20) 



A>0 



or, as the basis-independent expression 



p{r) = '^v^^^{r)u_^{r) - ]^vl{r)v^{r) 



duj 



P 



(4.21) 



where "tr" denotes a trace over Dirac matrices and P is a projection operator onto states orthogonal to the zero mode 
wave- functions. (We shall make use of the expression only where |r| — > oo and the zero-mode wave- functions have 
vanishing contribution. For this reason, explicit use of this projection will never be needed.) 

Case A5 = 0: Let us first consider the case where the axial vector gauge field is absent. In the asymptotic region, 
the Higgs field is 

if{r) = if,{r) + i^2(r) = cf, e'"« + ©(r'^), (4.22a) 
dMr) = -i^e^,^ </- e™*^ + 0(r-3). (4.22b) 

With 

VV^ = -df - dl + + a ■ d{ip, + ia^ip^), (4.23) 



we see that the derivatives of the Higgs field provide a long-ranged potential ^ 1/r in the Schrodinger equation (|4.12p . 
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It is easy to find the asymptotic behavior of the propagators in Eq. (|4.2ip by perturbative expansion in the deviation 
of the operator in Eq. (|4.23p from the free operator VD^ = —df — 5| + rri^. For example, 



1 



r = 



ip(j--r') I 



(4.24) 



2 



The right-hand side in this equation has support in the region where 0|r — r'| < 1 as it falls off exponentially with 
the distance |r — r'| when this distance is greater than 1/4). The second-term on the right-hand side of (|4.24p goes 
like 1/r and the corrections to it, represented by • • • , fall off faster and will not be needed. 
Using the asymptotic expression (|4.24p in Eq. (|4.2ip . we obtain 

p{r) = ^uL^(r)w_^(r) - ^v/,(r)z;,,(r) 

- 9 A I — ■ — ^ i^*id^ - idip*^) + ■■■ . 



Upon integrating this expression, we obtain 

Q = -i Index(ff) - ^y=^i^ / ' " i^V^V) (4-26) 

where, in the first term in the right-hand side, we have remembered that the number of zero modes is determined by 
the index and the line integral in the second term is taken on the circle at infinite radius. Using the index theorem (|4.9p 
and the asymptotic expression for ip in Eq. (|4.22ap . we obtain 

We have restored the possibility that /i could be negative in the first term by recalling that the sign of the energy of 
the zero modes is determined by fi. 

Case Ag ^ 0: The second case is when there is also an axial vector gauge field with asymptotic form 

Al{r)^-ne^^^ + 0{r-') (4.28) 

so that the covariant derivative of the Higgs field falls of at infinity 

{d + 2iAr^)ip{r) =0{r-^). (4.29) 

Here, we have assumed a power-law fall-off that is sufficiently fast for our purposes. In fact, for a classical field theory 
with a vortex solution, the covariant derivative falls off exponentially with distance from the vortex. Then, 



1 



PPt + ^2 + ^2 



J (27r)2 + 4^ + /i2 + cj2 



Corrections represented by • • • fall off at least as fast as l/r2 as r — > oo. The line integral in the phase factor is to be 
taken along a straight line between r and r'. (Since the axial magnetic field also goes to zero at least as fast as r~2 
as r oo, the path is not important for our purposes.) 

The trivial asymptotic form of (|4.30p means that the background fields do not contribute to the relevant asymptotic 
of the last term in (|4.2ip and the volume integral of that term vanishes. It therefore does not contribute to the total 
charge. We find that the charge in this case is entirely determined by the zero modes, 

1 n 
Q = -- Index(i/) sign(^) = sign(^). (4.31) 
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This is dramatically different from the result for the case without an axial vector gauge field quoted in Eq. (|4.27p . 
As we have seen, the difference can be attributed to the asymptotics of the background field configurations. Another 
way to understand it is to realize that, when the Higgs field approaches its asymptotic form its covariant derivatives 
as well as the axial magnetic field fall off quickly enough at r — ^ oo, so that stereographic projection can be used 
to map the problem of solving the Dirac equation on the plane to the problem of solving it on the sphere (where 
the vector field is a connection on a Wu-Yang monopole bundle). Then, the entire spectrum is discrete and, by the 
arguments following Eqs. (|4.14a[l and (|4.14bp wc can see that all non-zcro-mode solutions of the Schrodinger equation 
with eigenvalue A result in pairs of solutions of the Dirac equation: one positive E = -\/A + /i^ and one negative 
energy E = —^JX + jj? state. For this reason, only the zero modes can contribute to the spectral asymmetry and the 
contribution must be proportional to the index. 

As is well known, the axial gauge field is needed to render the vortex energy finite; it screens the infinite energy 
coming from the scalar field. Evidently, it also screens the irrational charge which arises from the staggered chemical 
potential. Some further insight into this phenomenon is given below. 

V. FIELD THEORY ANALYSIS 

An alternative method for finding the charge induced by the vortex background makes use of a field theoretic 
evaluation of the expectation value of the current in the "vacuum" state for the Dirac field operators in the vortex 
background, 

r{x) = {{i,{x)Yi^{x)) = -Tr [YS{x, x)] , (5.1) 

where S{x,y) is the Dirac field propagator for the Lagrange density (|2.3b[) . 

We consider first the theory without the axial gauge field and present Eq. (j2.3bp as 

£o=^(i7'^a,-$)V. (5.2a) 

where 

* = <Pi - i7V2 + 7V = <Pi - 17^2 + 7^3- (5-2b) 

Evidently, we need to invert 

S-\x, y) ^ -i {ij'-d, - $) S{x - y) . (5.3) 
This can be done pcrturbativcly in a gradient expansion for $. We set 

^{x) = M + S<i>{x), M = $(0), (5$(a;) = x''d^<^>{0). (5.4) 
In a graphical representation, a thick line denotes S{x, y) while a thin line represents the free propagator. 

^o(x,y)= /^-^e'^(^-), f^,^p^. (5.5) 



Hence we have 



6«I>(z) 
- i — 



h^{z) 80(z') 8a)(z) 80(z') 80(z") 

- N N + i H H ^ 

X z z' y X z z' z" y (5.6) 

with integration over the z-variables understood. The result of the calculation is 
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where 

3 

"^' = E^a = '/>' + M', "a = — , a =1,2,3. (5.7b) 



The second equality in Eq. (j5.7ap shows that J'^ is manifestly divergence-free. [The overall sign is set by the requirement 
that the isolated gap state is filled.] 

To evaluate the induced charge, we observe that the charge density for static fields is 

With our profile Eq. (|2.13p . this becomes 



Pl^) = 1 — 1 TT (5-9) 

47rr dr m{r) 



whose spatial integral yields 



= / d rp{r) = - 



2 m{r) 



(-isign(MO)) + i^K (5.10) 
2 ^ ' 2 m(oo) 



since the amplitude </) vanishes at the origin. This charge can be an irrational quantity, reducing to ±n/2 as the 
staggered chemical potential tends to ±0. 

It is noteworthy that the induced current (|5.7ap exhibits an 5*0(3) algebraic structure, even though neither the 
Lagrange density in Eq. (|5.2ap nor the propagator in Eq. ()5.3p put such structure into evidence. We shall explain 
below how this comes about. 

Another interesting point is that the current can take a simpler form after fields are redefined. First, we rewrite 
Eq. (j5.7ap in terms of ip and ip* 

J" - ^ e^"'' [/^ {d^^T [dpv) - ^^^Ji {ip*dp^ ~ ^dp^*)] . (5.11) 

Next, we define 

^ = 2mx^l-\x\\ M = "^(l-2|xP), (5.12) 
thereby expressing the current (|5.1ip as 

J. ^ Q^^. ^ i^'''"^^- - Xdpx*) ■ (5.13) 

This shows that J'^ is a total divergence and is manifestly conserved. 

Next we write the current when the axial gauge potential is present. Rather than calculating from first principles, 
we appeal to local axial gauge invariance, and promote all the derivatives in Eq. (jS.lip to covariant derivatives 

= d, + 21^5,. (5.14) 

However, the resulting expression is not conserved, but it can be made conserved by adding an axial gauge invariant 
term, which is linear in the axial gauge field. In this way we arrive at 

^'^ = a^''"'' {f^{D^vnDpV)-d^fi[^*{Di,^)~^{D^^*)]} + ^^F,-, (5.15) 

(Ji\ 1 1 L ^li fit 

where is the axial dual field strength 

F^^^le-^''F,^^=e'^^0d^A,f,. (5.16) 

As a check, the coefficient of the last term in Eq. (j5.15p can be computed from the relevant graph. When the axial 
gauge potential contribution to Eq. (|5.15p is separated, Eq. (|5.15p equals 

J" - 8^ {fiid^fridf^v) - d^ti[^*{df,v) - vidpv*)]} + [^A,^) . (5.17) 
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Therefore, the axial gauge potential's contribution to the charge density is 



which, for as in Eq. (j2.14p . equals 



and its contribution to the total charge is 



AQ= f d'r Apir) ^ ^ (5.20) 
J Z to(oo) 

since 05(00) = 1/2 while 05(0) = 0. This cancels the continuous dependence on ^(00) found in Eq. (|5.10p . leaving the 
same rational result obtained in the absence of the staggered chemical potential. 

Note that with variables defined as in Eq. (|5.12p . the current in the presence of the axial gauge field reads 



r 



or when the gauge field is separated 



= i^'^^^^o. [x*DpX - xDpX* - 21^5^] , (5.21) 



= i^''"''^" - - 2i (1 - 2|xn ^5^] . (5.22) 

Therefore, the total divergence feature and the conservation of the current are again explicitly exhibited also in the 
presence of the axial gauge field. In particular for the charge density, we have 

P = ^ ^''d,^ \x{D^x) - xiD^xf] + ^ e^'d^A,,. (5.23) 

Upon integration over space, the first term is cast on the circle at infinity, where the covariant derivatives of x vanish. 
The second term shows that the induced charge is exactly the vortex flux, equal to n/2 for as in Eq. (|2.14p . 

VI. INDUCED CHARGE FROM SYMMETRY ARGUMENTS 

As we observed in Sec.|Vl the form of the induced current (|5.7ap without an axial gauge field exhibits an 5*0(3) 
algebraic structure despite the absence of any such symmetry in the Lagrange density (|5.2aP i^ In this section we 
explain why this is so. Also we obtain expressions for the induced fractional charge from symmetry arguments. 

Our starting point is the Lagrange density (|2.3b[) with the mass terms collected into $, as in (|5.2bp . but written as 

C^=^ [^^ {id, + 75^5 J - mM.nJ ^ (6.1a) 

where 

Mi = a, M^^-if, A/3 = 73. (6.1b) 

The background fields mn^, functions on two-dimensional (2 + 1) dimensional space-time, are given by 

"1 = — , "2 = — 7 "3 = — ' (6-2) 

m m m 

which, since = -I- /i^, satisfy the local constraint 

l^nl + nl + nl^ri^. (6.3) 
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Despite the suggestive form in which the 3-dimensional vector n is written above, the Lagrange density (|6.1ap is not 
SU{2) symmetric because the matrices do not satisiy the su{2) algebra. 

However, it is the induced current (|5.7ap and not the starting Lagrange density (|6.1ap that exhibits the symmetry. 
Thus, let us turn our attention to the U{1) charge current induced by the background field mn, 

In Eq. (j6.4p . we are free to change integration variables as long as this transformation leaves the current unchanged. 
This we do through the nonunitary transformation 

fp = Xl5l^, "0 = X, (6-5) 
for some arbitrarily chosen constant unit vector TV. Thus, the induced current ()6.4p is now given by 



im,x] e'^''^^^ (xr-x) .... 

" f T^r- — 1 ifd^xc ^ 

/D[x,x] e'J 

where the transformed Lagrange density reads 

L^^X [r" {id, + ~ mE.iVJ X, (6.6b) 

with = 75 7"^ 7^ and = 75 7^ satisfying 

{r„r,,} = 2g,,„ [E,,E,] = 2ie,,,E„ [L,, EJ = 0. (6.6c) 
Since 75 does not commute with all E^, 

{Ei,7j = {E„7j = 0, {E3,75} = 21, [E3,75] = 0, (6.7) 
is an SU{2) singlet at = only 

A. Induced charge without axial flux 

In the absence of an axial gauge field, the Lagrange density (j6.6bp with = is an SU{2) singlet. The induced 
current and charges must therefore be 5*0(3) singlets given by 

r^C^ e-^f' e'^"^ d^n, 9^n„ (6.8a) 

I (fr J°(t,r)^C (6.8b) 
J 47r 

where Q, is the spherical angle (in units of An) covered by the mapping between the base space r S and a closed 
curve on the surface of the 2-sphere (j6.3p . to lowest order in a gradient expansion. Thus, we arrived at our previous 
result, except the factor C must still be determined. 

We can fix the constant C using our results for the fractional charge derived in the simple case when there is no 
staggered chemical potential: the charge is Q = — 1/2 when the midgap state is empty, and Q = 1/2 when the midgap 
state is filled. There is an ambiguity for the charge as given in Eq. (|6.8bp if the bound state is at exactly zero energy, 
because then it can be filled or empty, but this can be lifted by considering the case where ^ 0"^ with 0+ a positive 
infinitesimal. In this case, the bound state solution exists for an antivortex (n = —1), and E = — /i (see Refs. [20II21I 
and Sec. IIIip . so that the level is filled and the charge is thus Q = 1/2. 

Because ^ ^ O"*" or, equivalently, 0^, the spherical angle traced by the antivortex in the Ui 2 plane is just 

one full hemisphere (traced in the negative orientation): fl/Air = —1/2. Hence, the constant C = —1, leading to the 
induced current 

r^-j^ ^""^ e""' d^n, dpu, (6.9a) 



and the induced charge 



(6.9b) 
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B. Abelian formulation, including an axial flux 

Here we shall show that the induced charge in the presence of vortices in the ofF-diagonal masses and in an axial 
vector gauge potential is the same as that in a problem with constant off-diagonal mass and effective Abelian gauge 
flux. To this end, we make a further unitary transformation on the Lagrange density (|6.6bp 

X = X = C/^e (6.10a) 

The unitary matrix U is generated by the 4x4 matrices and is fixed by demanding that it takes the space-time 
dependent vector n in the fixed unit vector TV, 

(e.TV,) U{x) (S, n,{x)) U\x). (6.10b) 
It follows that U commutes with ^ but not with 75. With N = (0, 0, 1), this is achieved by choosing^ 

U{x) = e-'^ ^3 e+'^ ^2 e+'^ (6.11a) 
where the polar angle a and the azimuthal angle (3 parametrize n, 

n ~ (sin a cos (3, sin a sin /3, cos a) . (6.11b) 



The new Lagrange density reads 



where the matrix B,, is 



= e [F^ {id, + - 771S3] e (6.12a) 



B^ = U^^U^ A^, + md,U^. (6.12b) 

Now once the vector AT is fixed to (0, 0, 1), all the information about the original mass vortex and axial vector-gauge 
vortex is combined in B^^. The induced charge we want to compute is linear in these potentials (with higher orders 
suppressed by powers of m^^). Indeed to linear order, the current can only depend on the component of -B^ along 
the a = 3 direction. To sec this, consider a further rotation around the a = 3 direction by a constant angle 5, 

^ e+i|S3 5^e-'5S3. (6.13) 

The current is invariant under this rotation, but the components of B^, along the 0=1,2 directions do change. Hence, 
the induced current, at linear order, must not be a function of these components and it must depend solely on the 
component along the a = 3 direction 

K = Jtr (S3 

= itr [S3(C/75t/^A. + t/ia,C/^)] 

= ^9,/3-^(9,/3+2A5Jcosa. (6.14) 

We thus arrive at the result that the induced current and charge, computed using the Lagrangian ()6.12ap . are the 
same as those computed using the simpler Lagrange density 

A = Ur''(ia, + 6,I]3)-mE3]e (6.15) 

Finally, one last change of variables 

l^nY.^, e-?/, (6.16a) 

which does not affect the current, (again because of a trivial Jacobian in the path integral) and a redefinition of Dirac 
matrices 

% = S3F,, (6.16b) 
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which preserves their CHfFord algebra, gives the result that the induced current 

can be simply obtained from the Lagrange density with the gauge potential b^, and constant mass m 

= n[r i^d. + 15K) - H (6.17b) 

The flux due to is the only quantity left that retains any information on the mass and axial vortices, and thus 
it is the only variable controlling the value of the induced current and charge. The induced current must be an axial 
gauge invariant quantity and thus must be constructed from the axial flux due to b^^. The total charge, in particular, 
must be proportional to the total flux 

$5 = ^y d^r{dAb){r), (6.18) 

i.e., 

g = C$5. (6.19) 

The prefactor C is determined for the special case without a staggered chemical potential (// ^ 0+) and without an 
axial flux, in which case Q = +1/2. In this situation the polar angle a{r) — ^ 7r/2, cos Q;(r) 0, and b^{r) i d^(3{r) 
as r cx), so that for an antivortex the axial flux due to b^, is simply half the vorticity of the azimuthal angle f3(r): 
$5 = — ^. This flxes the constant C = — 1. We conclude that 

Q = -$5. (6.20) 

Equation (j6.20[) is the expression that we seek for the charge induced by static vortices in the mass and in the axial 
vector-gauge fields.— We now consider the following two cases (when /i > 0). 

(i) Static case with staggered chemical potential, no axial flux, and vorticity n in the mass, i.e., cos Q!(r) — > fj./m as 
r — > cxD and = 0: The flux is $5 = n ^ [1 — cos a(r — > 00)] and the induced charge equals 



= (6.21) 

2 V in/ 



(ii) Static case with staggered chemical potential and an axial vortex that screens the mass vortex, i.e., cosa(r) — > 
n/m and d^P{r) + 2 Agj^(r) -^0 as r 00: The last term in Eq. (|6.14p drops out and so does the dependence 
on the polar angle «(/•), along with the dependence on the staggered chemical potential. The flux due to b^^ is 
simply half the vorticity of the azimuthal angle P{r), $5 = ^ n, and thus the charge is pinned at the rational 
value 

Q = • (6.22) 



VII. SUMMARY 



The fractional charge induced by vortices supported by a complex- valued Higgs fleld carrying a [/(I) axial gauge 
charge of 2 that couples to massive Dirac fermions in (2 -I- l)-dimensional space-time was computed by three different 
techniques based on (i) the computation of a spectral asymmetry, (ii) a one-loop perturbative computation of the 
conserved fermion- number current, (iii) expressing the fractional charge in terms of an Abelian axial flux, respectively. 
The fractional charge can vary continuously as a function of an energy-reflection symmetry breaking parameter and 
thus can take irrational values. Remarkably, this fractional charge re-rationalizes to the value 1/2 taken in the presence 
of the spectral energy-reflection symmetry if an axial gauge field couples covariantly to both the Higgs and Dirac fields. 
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